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Abstract 



In Part I [25] of our study on 2D Euler equation, we established the spectral 
theorem for a linearized 2D Euler equation. We also computed the point spectrum 
through continued fractions, and identified the eigenvalues with nonzero real parts. 

In this Part II of our study, first we discuss the Lax pairs for both 2D and 3D Euler 
equations. The existence of Lax pairs suggests that the hyperbolic foliations of 2D and 
3D Euler equations may be degenerate, i.e., there exist homoclinic structures. Then 
we investigate the question on the degeneracy v.s. nondegeneracy of the hyperbolic 
foliations for Galerkin truncations of 2D Euler equation. In particular, for a Galerkin 
truncation, we have computed the explicit representation of the hyperbolic foliation 
which is of the degenerate case, i.e., figure-eight case. We also study the robustness 
of this degeneracy for a so-called dashed-line model through higher order Melnikov 
functions. The first order and second order Melnikov functions are all identically 
zero, which indicates that the degeneracy is relatively robust. The study in this 
paper serves a clue in searching for homoclinic structures for 2D Euler equation. The 
recent breakthrough result [24] of mine on the existence of a Lax pair for 2D Euler 
equation, strongly supports the possible existence of homoclinic structures for 2D 
Euler equation. 

PACS: 47.27 47.52 02.30.J 05.45 . 

MSG: 76 34 35 37 . 

Keywords: perverted heteroclinic orbit, higher order Melnikov function, degener- 
ate hyperbolic foliation. 
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1 Introduction 



In a series of works ([25], the current article, and [22]), we plan to build a dynamical 
system theory on 2D turbulence. The governing equation that we are interested in is 
the incompressible 2D Navier-Stokes equation under periodic boundary conditions. We 
are particularly interested in investigating the dynamics of 2D Navier-Stokes equation in 
the infinite Reynolds number limit and of 2D Euler equation. Our approach is different 
from many other studies on 2D Navier-Stokes equation in which one starts with Stokes 
equation to prove results on 2D Navier-Stokes equation for small Reynolds number. In 
our studies, we start with 2D Euler equation and view 2D Navier-Stokes equation for large 
Reynolds number as a (singular) perturbation of 2D Euler equation. 2D Euler equation is a 
Hamiltonian system with infinitely many Casimirs. To understand the nature of turbulence, 
we start with investigating the hyperbolic structure of 2D Euler equation. We are especially 
interested in investigating the degeneracy v.s. nondegeneracy nature of the hyperbolic 
foliation. Degeneracy means the coincidence of center-unstable and center-stable manifolds 
or unstable and stable foliations, i.e., figure-eight structure. If all the hyperbolic foliations 
of 2D Euler equation are of figure-eight, then there is no turbulent dynamics for 2D Euler 
equation. 

The recent breakthrough result [24] of mine on the existence of a Lax pair for 2D 
Euler equation makes my idea in the last paragraph even more realistic. The philosophical 
significance of the existence of a Lax pair for 2D Euler equation by the author, and more 
recently for 3D Euler equation by Steve Childress [4], is beyond the particular project 
undertaken here. If one defines integrability of an equation by the existence of a Lax pair, 
then both 2D and 3D Euler equations are integrable. More importantly, both 2D and 3D 
Navier-Stokes equations at high Reynolds numbers are near integrable systems. Such a 
point of view changes our old ideology on Euler and Navier-Stokes equations. 

In [25], we studied a linearized 2D Euler equation at a fixed point. The linear sys- 
tem decouples into infinitely many one- dimensional invariant subsystems. The essential 
spectrum of each invariant subsystem is a band of continuous spectrum on the imaginary 
axis. Only finitely many these invariant subsystems have point spectra. The point spec- 
tra can be computed through continued fractions. Examples show that there are indeed 
eigenvalues with positive and negative real parts. Thus, there is linear hyperbolicity. To 
understand the hyperbolic foliations of 2D Euler equation, especially the nature of degen- 
eracy v.s. nondegeneracy of the hyperbolic foliations, we first study a Galerkin truncation. 
The Galerkin truncation presented in this article allows us to compute the hyperbolic fo- 
liation explicitly. We call this Galerkin truncation the dashed-line model. The dashed-line 
model contains a control parameter e G [0, 1]. When the control parameter vanishes, the 
model decouples into a sequence of five- dimensional invariant subsystems. The eigenvalues 
of one of such invariant subsystems arc '"good" approximations of the eigenvalues of the 
original linearized 2D Euler equation. For this five-dimensional invariant subsystem, the 
hyperbolic foliation can be calculated explicitly with elegant formula representation. The 
hyperbolic foliation has the picture form of a "lip" . Both the upper and the lower halves 
of the "lip" are two-dimensional ellipsoidal surfaces. Orbits in each half of the "lip" are 
heteroclinic orbits spiraling into two fixed points as time approaches positive and negative 
infinities, and having a "turning point" ("perversion" [19] [13] [14] [15] [16] [18] [17]) in the 
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middle. We call these orbits perverted heteroclinic orbits. In terms of physical variables, 
such heteroclinic orbits represent heteroclinic profile-evolutions approaching two different 
stationary profiles as time approaches positive and negative infinities. Recently, there are 
more and more interests in studying heteroclinic and homoclinic orbits in fluids [8] [6] [7]. 
For the two-dimensional Kelvin- Helmholtz problem, Craig and Groves found a homoclinic 
orbit approaching a temporally periodic profile as time approaches positive and negative 
infinities for a fourth order normal form [7]. Numerical calculation indicates that the eigen- 
values with zero real parts of the dashed-line model when the control parameter vanishes, 
form four narrow bands. When the control parameter is increased from to 1, the sizes 
of these four bands increase and congregate into the one band continuous spectrum of the 
original linearized 2D Euler equation. In this sense, the eigenvalues with zero real parts 
of the dashed-line model approximate the continuous spectrum of the original linearized 
2D Euler equation. We are interested in investigating the degeneracy v.s. nondegeneracy 
nature of the hyperbolic foliation for 2D Euler equation. As a first step, we will study the 
persistence of the "lip" structure for the dashed-hne model through Melnikov functions. It 
turns out that both the first order and the second order Melnikov functions are identically 
zero. 

There have been a lot of works on Melnikov functions since the original work of V. K. 
Melnikov [27]. It has been proved rigorously that a Melnikov function is the leading order 
term of a certain signed distance between center-unstable and center-stable manifolds, see 
for example [20] [33] [5] [26]. Typically, a Melnikov function is a temporal integral from 
negative infinity to positive infinity of the inner product of the gradient of certain invariant 
with the perturbation term evaluated along an unperturbed heteroclinic or homoclinic orbit, 
and this integrand often has exponential decay property as time approaches positive and 
negative infinities. The importance of Melnikov functions is to determine the intersection 
between the center-unstable and the center-stable manifolds. The argument is as follows. 
Often the Melnikov function is easily computable. First one seeks zeros for the Melnikov 
functions; then implicit function theorem implies nearby zeros for the signed distance. 
The intersection between the center-unstable and the center-stable manifolds implies the 
existence of orbits homoclinic to the center manifold. Other variety of Melnikov functions 
include subharmonic Melnikov functions [20], which are temporal integrals over the finite 
period intervals, and exponentially small Melnikov functions [12], etc. In some special 
circumstance, the Melnikov function can be identically zero as a function of parameters. In 
such case, one needs to calculate the next leading order term of the signed distance, called 
the second order Melnikov function. If the (n — l)-th order Melnikov function is identically 
zero as a function of the parameters, the n-th order Melnikov function serves as the leading 
order term of the signed distance. There have been some works on second order Melnikov 
functions, see for example [29] [30] [34] [28] [9] [21]. In this paper, we are going to derive 
high order Melnikov functions without rigorous justification for the dashed-line model, and 
calculate them to study the persistence of the figure-eight structure. 

Our purpose of introducing the dashed-line model is for investigating the hyperbolic 
foliation of 2D Euler equation, which is different from those of other models. For example, 
the shell model [3] [32] models the interaction of energy in different shells in the spectral 
space. The difference between the dashed-line model and the shell model can be summarized 
as follows: 1. The nonlinear terms in the dashed-hne model are derived from Galerkin 
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truncations, while the nonhnear terms in the shell model are postulated. On the other hand, 
both the dashed-line model and the shell model conserve the energy and the enstrophy. 2. 
The dashed-line model captures the eigenvalues of the linearized 2D Euler equation, while 
the shell model has no such information at all. 3. The dashed-line model is claimed to 
model the dynamics of the 2D Euler equation in a neighborhood of a line of fixed points and 
hopefully also the global dynamics. The variables in the dashed-line model have precise 
physical meanings. The shell model is claimed to model the transfer of energy in different 
spectral shells. The variables in the shell model are artificial. 

The article is organized as follows: In Section 2, we discuss the Lax pairs for both 
2D and 3D Euler equations. In Section 3, we will introduce the dashed-line model and 
discuss its rationality. In Section 4, we will compute explicitly the figure-eight structure 
for the dashed-line model at the special parameter value e = 0. In Section 5, we study 
the persistence of the figure-eight structure when e 7^ through higher order Melnikov 
functions. Section 6 is the conclusion. 

2 Lax Pairs for 2D and 3D Euler Equations 

In this section, we are going to discuss the most recent results on the Lax pairs of 2D and 
3D Euler equations. 

2.1 The Lax Pair for 2D Euler Equation 

We consider the two-dimensional incompressible Euler equation written in vorticity form. 



dQ _ dQ dQ 
dt dx dy 

(2.1) 

du dv ^ 
dx dy ' 

where Vt is vorticity, u and v are respectively velocity components along x and y directions. 
If we define the stream function ^ by, 

_ _d^ 
dy ^ dx ^ 

then we have the relation between vorticity Vt and stream function ^, 

^ dv du 

Q = = . 

dx dy 

The 2D Euler equation (2.1) can be written in the equivalent form, 

^ + {*,Q} = 0, (2.2) 
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where the bracket { } is defined as 




{f,9} = {dj){dyg)-{dyf){d,g) . 
Theorem 2.1 ([24]) The Lax pair of the 2D Euler equation (2.2) is given as 

,0, 

where 

L^ = {n,ip}, Aip = {^,ip}, 

and X is a complex constant, and ip is a complex-valued function. 

The compatibihty condition of the Lax pair (2.3) gives the 2D Euler equation (2.2), i.e. 

dtL = [L,A] , 

where [L,A] = LA — AL, gives the Lax representation of the 2D Euler equation (2.2). In 
investigating 2D Euler equation through the Lax pair, two areas seem promising. The first 
area is Darboux transformation. Darboux transformations have been utilized for generating 
explicit representations for homoclinic structures for soliton equations, see for example [23]. 
Up to now, the Darboux transformation for the above Lax pair has not been found. On the 
other hand, in the present paper, we do not use the approach of Darboux transformation for 
investigating homoclinic structures for 2D Euler equation, rather use the classical approach 
of Galerkin truncations. Investigating the Darboux transformation is our future research. 
The other area is inverse scattering transform. Inverse scattering transforms have been 
exploited for solving the Cauchy problems of soliton equations, see for example [11] [2]. In 
building the inverse scattering, it is crucial to have constant coefficient differential operators 
in (especially the spatial part of) the Lax pair. For that reason, one may start with the 
following Lax pair [35] , 

/ Diip + {Q,(p} = X(p , , . 

\ dtip + D2ip + {S,ip}^0 , ^^-^^ 

where 

7^ 9 n d ^ d 

Di = a— + /3—, 1)2 = 7— + (5—, 
ox oy ox oy 

a,f3,'y, and 6 are real constants, A is a complex constant, S* is a real-valued function, and 
(y9 is a complex-valued function. The compatibility condition of this Lax pair gives the 
following equation instead of the 2D Euler equation, 

— + {S, n} + L>2Q - DiS = 0. (2.5) 

A promising approach is that one can first build the inverse scattering for (2.4) and (2.5), 
and then take the limits 
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to get results for 2D Euler equation. A special case of the Equation (2.5) is the following 
system of equations [35], 

' § + {s,n} = Q, 

< 

DiS = . 

In this paper, we are not going to follow the approach of inverse scattering transform 
either. Investigating the inverse scattering transform is our future research. Here the Lax 
pair structures make it plausible for the existence of homoclinic structures for the 2D Euler 
equation. 

2.2 The Lax Pair for 3D Euler Equation 

We consider the three-dimensional incompressible Euler equation written in vorticity form, 

dtVl + {u ■ V)fi - (fi ■ V)m = , (2.6) 

where u = {ui,U2,Us) is the velocity, Q = {Qi,Q2,^3) is the vorticity, Q = V x u, and 
V ■ M = 0. u can be represented by Q for example through Biot-Savart law. More recently, 
to the author's surprise, the Lax pair for 3D Euler equation was established by Steve 
Childress. 

Theorem 2.2 ([4]) The Lax pair of the 3D Euler equation (2.6) is given as 
where 

L(p — Q, ■ \/(p — (p ■ VJl , A(p — u • V</7 — (/? • , 
\ is a complex constant, and Lp = {ipi,ip2,(p3) is a complex 3-vector valued function. 

The compatibility condition of the Lax pair (2.7) gives the 3D Euler equation (2.6), i.e. 

dtL^[L,A] , 

where [L, A\ = LA — AL. gives the Lax representation of the 3D Euler equation (2.6). 

Unfortunately, when one does 2D reduction to the Lax pair (2.7), one gets L = 0. 
Therefore, the Lax pair (2.7) does not imply any Lax pair for 2D Euler equation. 

In investigating 3D Euler equation through the Lax pair, especially the question on the 
possibility of finite time blow up solutions, the two areas: (1). Darboux transformations, 
(2). inverse scattering transforms, also seem promising. Up to now, the Darboux transfor- 
mation for the Lax pair (2.7) has not been found. Investigating the Darboux transformation 
and the inverse scattering transform for the Lax pair (2.7) is our future research. 

In building the inverse scattering, it is crucial to have constant coeficient differential 
operators in (especially the spatial part of) the Lax pair. For 3D Euler equation, we may 
start with the following Lax pair 
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Proposition 1 We consider the following Lax pair, 



j Lip = Xip , , . 

where 

Lip = Q ■ V(/9 — ■ Vr2 + Diip , A(p = q ■ V<f — <f ■ Vg + D2ip , 

X is a complex constant, ip = {pi,p2,'^3) is a complex 3-vector valued function, q = 
(gi, q2, qs) is a real 3-vector valued function, Dj = a^^^ ■ V, (j = 1,2), a^^^ = {a[^\ \ a^^) 
are real constant 3-vectors. The compatibility condition of this Lax pair gives the following 
equation instead of the 3D Euler equation, 

dtVt + {q ■ V)n - {n ■ V)g + - D^q ^ . (2.9) 

A specialization of (2.9) is the following system of equations, 

' dt^ + {q ■ v)n -{n-v)q = o , 

< 

_ Diq = L>2Q . 

Proof: The proof of this proposition is a trivial direct verification. □ 
A promising approach is that one can first build the inverse scattering for (2.8), and 
then take the limits 

Q,W^O, (j = l,2) , q^u, 

to get results for 3D Euler equation. Here the Lax pair structures make it plausible for 
the existence of homo clinic structures for the 3D Euler equation based upon the study on 
soliton equations [23]. 

Remcirk 2.1 // one defines integrability of an equation by the existence of a Lax pair, 
then both 2D and 3D Euler equations are integrable. More importantly, both 2D and 3D 
Navier-Stokes equations at high Reynolds numbers are near integrable systems. 

3 Introduction of the Dashed-Line Model 

We consider the two-dimensional incompressible Euler equation written in vorticity form(2.1)| 
under periodic boundary conditions in both x and y directions with period 27r. We also 
require that both u and v have means zero, 

/'27r /'27r rliv rliv 

u dxdy — j V dxdy — 0. 
Jq Jo Jo Jo 

We expand fl into Fourier series, 

feez2/{o} 
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where uj_^. = uJf: , k = {ki, k2)'^, X = {x,y)'^. In this paper, we confuse with (0,0)'^, the 
context will always make it clear. By the relation between vorticity Q and stream function 
^, the system (2.1) can be rewritten as the following kinetic system, 



1 ' 



where A{p, q) is given by. 



(3.1) 



Mp^Q) = 2^\q\ ^-\p\ ^]{Pi<l2 - P2qi) 



1-2 



\P\ 



Pi Qi 

P2 q2 



(3.2) 



where |gp = 9i + for q = {qi, q2)^, similarly for p. 

For any two functionals Fi and F2 of {cok}, we define their Lie-Poisson bracket as 



k+p+q=0 



Ql Pi 
q2 P2 



dFi dF2 
dUJ^ dUJ^ 



(3.3) 



Then the 2D Euler equation (3.1) is a Hamiltonian system [1], 



(3.4) 



where the Hamiltonian H is the kinetic energy, 

1 



^ fcez2/{o} 



(3.5) 



Following are Casimirs (i.e. invariants that Poisson commute with any functional) of the 
Hamiltonian system (3.4): 



Jn — E '^'=1 ■ ■ ■ ^kn- 

kl + -+kn=0 



(3.6) 



3.1 Preliminaries on Linearized 2D Euler Equation 



In this subsection, we discuss the preliminary results on linearized 2D Euler equation known 
from [25]. 

We denote {oJk}kez^/{o} by OJ- We consider the simple fixed point uj*: 



r, = 0, if A; 7^ p or - p. 



(3.7) 
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of the 2D Euler equation (3.1), where F is an arbitrary complex constant. The linearized 
two-dimensional Euler equation at lu* is given by, 

ujk = A{p, k-p)r Uk-p + A{-p, A; + p) f Uk+p ■ (3.8) 
Definition 1 (Classes) For any k G Z^/{0}, we define the class to be the subset of 

zy{o}: 

n & Z, p is specified in (3.7) >. 

Sec Fig. 3.1 for an illustration of the classes. According to the classification defined in Defi- 
nition 1, the linearized two-dimensional Euler equation (3.8) decouples into infinitely many 

invariant subsystems: 

^k+np = Mp, k + {n- l)p) F 

(3.9) 

+ A{-p, k + {n + l)p) F (^fe+(„+i)p . 



Theorem 3.1 The eigenvalues of the linear operator Cf, defined by the right hand side of 
(3.9), are of four types: real pairs (c,—c), purely imaginary pairs (id,— id), quadruples 
(±c±id), and zero eigenvalues. 

The eigenvalues can be computed through continued fractions. 

Definition 2 (The Disk) The disk of radius \p\ in Z^/ {0}, denoted by D^p^, is defined as 

L'l,l = |A;eZ7{0} I 1^1 < 
The closure of D\p\, denoted by D\p\, is defined as 

D\p\^i^k&Z^/{Q} |A;|<|p||. 
See Fig. 3.1 for an illustration. 

Theorem 3.2 (The Spectral Theorem) We have the following claims on the spectra of 
the linear operator Cf.: 

1. IfT^f.P[D\p\ = 0, then the entire £2 spectrum of the linear operator Cj^ is its continuous 



spectrum. See Figure 3.2, where b = —^\T\\p\~ 



pi ki 
P2 k2 



2. If n Z)|p| 7^ 0, then the entire essential £2 spectrum of the linear operator Cf. is its 
continuous spectrum. That is, the residual spectrum of C^. is empty, crj.(>C^) = 0. The 
point spectrum of C^, is symmetric with respect to both real and imaginary axes. See 
Figure 3.3. 
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3.2 Rationality of the Dashed-Line Model 

To simplify our study, we study only the case when cuk is real, VA; e Z^/{0}, i.e. we only 
study the cosine transform of the vorticity, 

Q — LOk cos(/c • X) , 

and the 2D Euler equation (2.1;3.1) preserves the cosine transform. To further simplify 

our study, we will study a concrete dashcd-linc model based upon the line of fixed points 
(3.7) with the mode p = (1, 1)"^ parametrized by F. When F 7^ 0, each fixed point has 4 
eigenvalues which form a quadruple. These four eigenvalues appear in the only unstable 
invariant hnear subsystem labeled by A; = (—3, —2)'^. We computed the eigenvalues through 
continued fractions, one of them is [25]: 

A = 2A/|r| = 0.24822302478255 + i 0.35172076526520 . (3.10) 

Sec Figure 3.4 for an illustration. We hope that a Galerkin truncation with a small number 
of modes including those inside the disk D\p\ can capture the eigenvalues. We propose the 
Galerkin truncation to the linear system (3.9) with the four modes k + p, k + 2p, k + 3p, 
and k + Ap, 
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i2\b\ 




* 




-i2\b\ 



K{A} 



Figure 3.2: The spectrum of £^ in case 1. 



u;2 = AiFo;! - A3ra;3 , 
Us = A2ra;2 - A4ra;4 , 

Prom now on, the abbreviated notations, 

'^n = cc;^+„p, An = A{p,k + np) , Am,n = A{k + mp,k + np) , (3.11) 

will be used. The eigenvalues of this four dimensional system can be easily calculated. It 
turns out that this system has a quadruple of eigenvalues: 

r 



A = ±- 



= ± 



2V10 
T 



(^-^ X 0.7746 X e^'^' , 



(3.12) 



where 9i = arctan(0.845), in comparison with the quadruple of eigenvalues (3.10), where 

A= ± (^^^ X 0.43 X e±*^2 , 

and 62 = arctan(1.418). Thus, the quadruple of eigenvalues of the original system is 
recovered by the four-mode truncation. Along the above line of thinking, we "chop" the line 
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i2\b\ 



-i2\b\ 



K{A} 



Figure 3.3: The spectrum of Ct in case 2. 



T,^ as follows, 



= -A 



5j+2 



r a;5j- 



+2 , 



= A 



5j+3 



r UJ^j 

r u;5j- 



;+2 



+3 



That is, we "chop" the hue T,^ at the points n — 5j, Vj G Z. When j ^ 0, this system has 
two complex conjugate pairs of purely imaginary eigenvalues. 



A = ±i(r/V2)i 



where 



6 = 



-^5j+1^5j+2 



^5j+2^5j+3 



± 



C — ^5^-1.1^5^+2^5^-1-3^5^+4 , 

6<0, OO, 0<6^-4c<6^. 



When j — > ±oo. 



A ^ ±i(r/2)i 



/3 + 5 V5 



± 
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K{A} 



Figure 3.4: The quadruple of eigenvalues of the invariant system labeled by ^ = (—3, —2)^ 
when p — (1, 1)^. 



= ± i 0.31 r , ± i 0.8 r . 

For example, when j — 1, 

X = ± i 0.2937609 T , ±i 0.7736967 T ; 

when j — 2, 

A = ± i 0.3057701 T , ±i 0.8007493 F . 

For all j G Z, the distribution of the eigenvalues is illustrated in Fig. 3. 5. To recover the 
original linear system from the "chopped" system, we introduce the homotopy parameter 
e G [0, 1], and consider the homotopy system, 

where 

^ ^ r 1, iin^5j,\fjeZ , 
^ 1 e , if n = 5j , for some j E Z . 

The homotopy deformation of the spectra for this homotopy system is calculated numeri- 
cally by Thomas Witelski and shown in Fig. 3.6 for a special value of F. As e is increased 
from to 1, the quadruple of eigenvalues moves along a curve, and the sizes of the four 
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Figure 3.5: The distribution of tlie eigenvalues of tlie cliopped system. 



Figure 3.6: The homotopy deformation of the spectra for the homotopy system. 



bands of eigenvalues with zero real parts increase, and finally these four bands congregate 
into one band when e = 1 which is the continuum spectrum. From such numerical calcula- 
tion, one can clearly see the homotopy transition of the spectra from a decoupled "chopped" 
subsystem to the original linearized 2D Euler equation. From the above discussions, it is 
natural to propose the following dashed-line model , 

(3.13) 



Up 



n&Z 



to model the hyperbolic structure of the 2D Euler equation, connected to the line of fixed 
points (3.7) with p = (1, 1)"^. Figure 3.7 illustrates the collocation of the modes in this 
model, which has the "dashed-line" nature leading to the name of the model. 

The dashed-hne model has the same properties as 2D Euler equation, of conserving the 
kinetic energy and enstrophy, and being a Hamiltonian system with the same Lie-Poisson 
bracket structure. For any two functionals Fi and F2, we define their Lie-Poisson bracket 
as follows. 



k+q+r=0 



ri qi 
r2 ?2 



dFi dF2 

^k^kT\ 7^ 



The kinetic energy is given by. 



nez 
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The dashed-line model (3.13) is a Hamiltonian system with the kinetic energy as the Hamil- 
tonian, 

< 

, Up = {up, H} . 

The enstrophy 

is still a constant of motion. The invariance of kinetic energy and enstrophy is a consequence 
of the simple relations, 

jif2A{p, q) + k) + = , 

A{p,q) + A{q,k) + A{k,p) = Q ■ 

whenever k + p + q = 

The dashed-line model is also a homotopy system parametrized by the homotopy pa- 
rameter e e [0,1]. When e = 0, the figure-eight hyperbohc structure of the dashed-hne 
model is computable as studied in next section. As a first step toward understanding the 
homotopy deformation of such hyperbolic structure, we are going to use higher order Mel- 
nikov functions to study the persistence of such figure-eight hyperbolic structure when e is 
small. 
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4 The Degenerate Hyperbolic Foliations of the Dashed-] 
Line Model When e = 

When e = 0, the dashed-Une model takes the form: 

<^5j+3 = ^5j+2<^p<^5j+2 — A5j+4UpL05j+4 , {j G Z) , (4-1) 



P 



-E 



^5i+l,5i+2t^5i+l'-^5i+2 + A5j+2,5j+3^5j+2i^5j+3 



+ ^5j+3,5j+4<^5j+3<^5j+4 

and the equation for the decoupled variable cu^j is given by, 

This system has a sequence of invariant subsystems, one of which carries the hyperbohc 
structure. 

4.1 Invariant Subsystems When e = 

For each fixed j e Z, we have the seven dimensional invariant subsystems: 

^5j+l — — ^5i+2'^p<^5i+2 , 

f^Sj+a = ^5jr+2^pl^5j+2 — A^j^4^U!pU!5j+4^ , (4-2) 

= — [^5i+l,5j+2'^5i+l'^5i+2 + ^5i+2,5i+3'^5i+2'^5i+3 
+ ^5i+3,5i+4<^5i+3'-^5i+4] , 

and the equations for the decoupled variables cu^j and a;5j+5 are given by, 

(jj^j — —A^j^iOjpOj^j+i , 

^5j+5 — ^5j+4<^p<^5j+4 ■ 



The kinetic energy 



1 



1 



E-"- 2 I 2 



and the enstrophy 



i<;<4 

are constants of motion for this system. 
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4.2 Perverted Heteroclinic Orbits When e = 



Prom the studies in last subsection, the {j = 0) hnearized system at the fixed point (3.7) 
has a quadruple of eigenvalues with nonzero real parts, and all {j ^ 0) linearized system at 
the fixed point (3.7) has two pairs of purely imaginary eigenvalues. Next we focus on the 
invariant subsystem (4.2) when j = 0, and generate explicit expressions for the unstable 
and stable manifolds of the line of fixed points (3.7). Notice that 



3 1 



4 

^1,2 = Ai-A2^--, A2,3 = , ^,4 = -^1,2 ; 

o 



(4.3) 



then the invariant subsystem (4.2) for j = can be rewritten as follows: 

COl = —A2 Up OJ2 , 

(02 = Ax UJp uJi - A2 UJp uz , 

ujz = A2 UJp UJ2 — Ai UJp uji , (4.4) 

a;4 = A2 UJp UJ3 , 

OJp = ^1,2 (<^3 <^4 - UJl UJ2) , 

and the equations for the decoupled variables ujq and uj^ are given by, 

uo = -Ai UJp UJl , 
UJ5 = Ai UJp UJ4 . 

There are three invariants for the system (4.4): 

I = 2Ai^2iuiUJ3 + UJ2UJ4) + A2UI , (4.5) 
U = Aiiujf + ujI) + A2{ujI + ujD , (4.6) 
J = ujI + ujI + uj^ + ujI + ujI . (4.7) 

J is the enstrophy, and U is a linear combination of the kinetic energy and the enstrophy. 
/ is an extra invariant which is peculiar to this invariant subsystem. With /, the explicit 
formula for the hyperbolic structure can be computed. 

The common level set of these three invariants which is connected to the fixed point 
(3.7) determines the stable and unstable manifolds of the fixed point and its negative —uj*: 

ujp^-r,ujn^O (neZ). (4.8) 

Using the polar coordinates: 

UJl — r cos 9 , UJ4 — r sin 9; uj2 — p cos , UJ3 — p sin 1? ; 
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we have the following explicit expressions for the stable and unstable manifolds of the fixed 
point (3.7) and its negative (4.8) represented through perverted heteroclinic orbits: 



CO, 



— r tanhr, 



A, 



A2-A1 

A2 



r sech T , 



2k 



In cosh T + 9o , 



(4.9) 



9 + '& 



-Ai 

— arc sin 
TT + arc sin 



1 rjz 

2\l -Ai 



1 / M 

2 V -Ai 



(«: > 0) , 
(«: < 0) , 



where Ai and A2 are given in (4.3), r = /tFt+ro, (tq, ^^o) are the two parameters parametriz- 
ing the two-dimensional stable (unstable) manifold, and 



K 



-A^A2 cos{e + ^) = ±J-A,A2\ 1 + 



4^1 



The two auxilliary variables cuq and cu^ have the expressions: 



U}5 = 



sech T |sin[/3 In cosh t -\- 9o\ — — cos[/3 In cosh r -|- ^0] | , 
sech T I cos[/3 In cosh t -|- ^0] + 4 sin[/3 In cosh r -|- ^0] I , 



/3 



where 



Remark 4.1 /n this remark, we will discuss the complex version of the system (4-4) '■ 



OJl = -A2 uJ^ OJ2 , 

UJ2 = Ai Up LOi — A2 LOs , 

Us = A2 Up U2- AiU^ U4, 

U4 = A2 Up U3 , 



(4.10) 



Ur, 



Ai^2 (<^3 <^4 - ^i^l UJ2) 



This is a ten dimensional Hamiltonian system with no reality restriction. For this system 
(4-iO), we have one complex invariant 



Ic = 2Ai^2{uJlU3 + U2U4) + A2UI , 
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and two real invariants 

U, = Ai{\uj^\^ + \uj4^) + A2{\UJ2\^ + \ujs\^) , 
Jc = |<^pP + |<^l|^ + |<^2|^ + |<^3|^ + ^4!^ ■ 

Thus, we have total four invariants which are not enough for the integrability of the Hamil- 
tonian system (4-10). When we consider the reality restriction, we have three invariants. 
Their common level sets in the five dimensional real phase space, are two-dimensional. The 
stable and unstable manifolds of the fixed point are also two-dimensional, and these three 
invariants are enough to determine the stable and unstable manifolds. 

4.3 Discussion on the Perverted Heteroclinic Orbits 

The oscillatory nature of the heteroclinic orbit reflects the fact that the quadruple of eigen- 
values with nonzero real parts, also has nonzero imaginary parts. The interesting nature of 
the heterochnic orbits is that they all have a "perversion" which refers to the neighborhood 
of the turning point. The portions of the heteroclinic orbit in the neighborhoods of u* and 
—uo* can be viewed as oriented helices. Then the entire heteroclinic orbit is a connection 
between a right-handed helix and a left-handed helix, and the connection part has to be a 
perversion. See Figure 4.1 for an illustration. Prom the study on the hnearized 2D Euler 
equation [25], we realize that the eigenvalues of the fixed point uj* only depend upon the 
modulus \ujp\ (the same fact is true for the dashed-line model too). Therefore, uj* and —uj* 
have the same eigenvalues. In fact, the eigenvalues appear in quadruples which make the 
connection of a right-handed helix with a left-handed helix feasible. 

The term "perversion" was given by the 19th century topologist Listing to describe the 
spontaneous switching of a helical structure of one handedness to its mirror image. Tendril 
perversion in climbing plants had been described at length by Charles Darwin in his book 
[10]. Darwin interpreted the tendril perversion using the language of elasticity. Complete 
mathematical study on elastic filaments modelling tendril perversion has been developed 
by A. Goriely and M. Tabor [19] [13] [14] [15] [16] [18] [17]. Figure 4.2 shows a climbing 
plant with one tendril perversion [31]. One tendril can have many perversions [10], see 
Figure 4.3. Perversions can occur in problems like the microscopic properties of biological 
fibers such as cotton or the formation of bacterial macrofibers, telephone cords (Figure 
4.2), false-twist technique in the textile industry etc. [19]. The explicit expression (4.9) 
of the perverted heterochnic orbits shows that the unstable manifold 14^" (cu*) of the fixed 
point UJ* (3.7) is the same as the stable manifold W^{uj*) of the negative —uj* (4.8) of uj*, 
and the stable manifold of the fixed point uj* (3.7) is the same as the unstable manifold of 
the negative —uj* (4.8) of uj*. Both W'^{uj*) and W^{uj*) have the shapes of "painted eggs" 
(Figure 4.4). W^{uj*) and W^{uj*) together form the "lip" (Figure 4.4) which is a higher 
dimensional generalization of the heteroclinic connection on plane. As the first step toward 
understanding the degeneracy v.s. nondegeneracy nature of the hyperbolic structure of 
2D EuIcr equation, we are interested in the e-homotopy deformation of such hyperbolic 
structure for the dashed-line model. In the next section, we will use higher order Melnikov 
functions to study the "breaking" or "persistence" of such hyperbolic structure when e is 
small. 
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Figure 4.1: The connection part of a right-handed heUx and a left-handed helix has to be 
a perversion. 



Figure 4.2: (a). Tendril perversion in Bryonia dioica. Illustration from Sachs' Text-book of 
Botany (1875); (b). Perversion in a telephone cord (Goriely and Tabor, 1997). The effect 
is achieved by fully stretching out and untwisting the cord and slowly bringing the ends 
together. The above figures are taken from Goriely and Tabor's papers (see References). 



Figure 4.3: A caught tendril of Bryonia dioica, spirally contracted in reversed directions. 
Illustration from Darwin's book: The Movements and Habits of Chmbing Plants (1875) 
(see References). 
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Figure 4.4: The unstable and stable manifolds of the fixed point ui* (3.7): (a) and (b) show 
two "painted eggs" and (c) shows a "lip" . 

5 The Higher Order Melnikov Functions 

In this section, wc are going to study the persistence of the hyperbolic structures given 
by (4.9), i.e., the persistence of the "painted eggs" and "lips" as shown in Figure 4.4 of 
the invariant subsystem (4.4) in the dashed-line model (3.13) when e 7^ 0. We will use 
Melnikov functions to detect such persistence, and the Melnikov functions will be built 
upon the following two invariants of the dashed-line model (3.13) when e = 0, which are 
actually invariants of the invariant subsystem (4.4): 



where F is a linear combination of the invariants / and J given in (4.5) and (4.7), V — 



From the spectral distribution Figure 3.5, when e = 0, the fixed point (3.7) or (4.8) has 
codimension 2 center- unstable and center-stable manifolds W'^'^{-^T) and iy^*(±r) under 
the flow (4.1), and the "painted eggs" are two-dimensional submanifolds of them. For 
sufficiently small e, W^™(±r) and W'i^T) perturb into W^''{±T) and W^f (iP). We are 
going to use the new notation 



U = Ai{luI + Lul) + A2{CU2 + i^l) , 

V = A2{ujI + + ul + ul) 



(5.1) 



(5.2) 



A2J - 1. 




(5.3) 



(5.4) 



(5.5) 



5.1 The Melnikov Functions 



The (first order) Melnikov functions are given by: 




(5.6) 
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J — oo 



My = 



/oo 
{^V,g)lio) dt 
-oo 



J —oo 



/ [2AoA2U)pOJoUJi — 2A0A2UJPOJ4UJ5 



-2AoAi^2COpUJoUJ3 + 2^o^l,2<^p<^2<^5] L(o) dt 



- 9A r (O).(O) (0).(0)^ 

J— 00 lAi 



since the second round bracket is an odd function in r. The above representations will be 
reproduced in the process of deriving higher order Melnikov functions in the next subsection. 
The above calculations show that the Melnikov functions are identically zero as functions 
of and r. Thus the separation between W^"'{-^T) and W^'^{-^T) is at least of order 0{e^). 
To further detect the separation between W^'^{±V) and ^^^'^'^(ir), we need to compute the 
second order Melnikov functions which are the leading order terms of the signed separation 
distances. 

5.2 The Derivation of Higher Order Melnikov Functions 

The derivation in this subsection is given without rigorous justifications, and rigorous 
justifications are future works. 

Let ijj^^\t) denote a heteroclinic orbit given by (4.9), and let S be a codimension 1 
hypersurface which is transversal to the orbit cj'^°'(t) at uj^^\{)). See Figure 5.1 on the setup 
of this derivation. Without loss of generality, we assume that in backward time a;*^°)(t) 
approaches to* (3.7) and in forward time uj^'^\t) approaches —u* (4.8). The submanifolds 
iycs^_P^ n E and W^^{—T) fl E have codimension 2 in E, and the codimensions are co- 
ordinated by VU and VV^. Let uj^'^\t,e) denote an orbit of the (e 7^ 0) system (3.13) 
which approaches uj* in backward time. Let C(j*^")(0, e) be the intersection point of cu'^'^^t, e) 
with E. On the submanifold W^^{—T) n E, let a;^*''(0, e) be the point which has the same 
coordinates with a;*^"^(0, e) except the Vt/ and W directional coordinates in E. Denote by 
u;^^\t, e) the orbit of the (e 7^ 0) system (3.13) with initial point u;^^\0, e). 

Define the signed distances: 



du = (Vf/(0),cuH(0,e)-c.;W(0,e)) 

= (Vt/(0),cu(")(0,e) -cu(°)(0)) 

dy = (Vy(0), a;(")(0,e) -cu(^)(0,e)) 

= (Vy(0),a;("'(0,e) -u;(°)(0)) 



(5.7) 
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Figure 5.1: The geometrical setup for deriving the second order Melnikov functions. 



-(V\/(0) , J'\0, e) - J^\0)) . (5.8) 

From now on, we will do the derivation for du, and the derivation for dy is the same. We 
define the following two functions: 

AUt,e) = (VC/(i),a;(^)(i,e)-a;(°)(0), te[0,oc), 
A^it,e) = {WU{t,e),J-\t,e)-J'\t)), ie(-oo,0], 

where VU{t) = VU{uj^^\t)). Then 

= A^(0,e)-A+(0,e). (5.9) 

Next we will do the derivation for A^(t, e) and the derivation for A^(t, e) is the same. For 
any T > 0, A^(i, e) as a function of e has the Taylor expansion: 

oo 

A^(t,e) = ^e"Af;'"H^) , te[0,T], (5.10) 

n=l 

where 

and uj^'\t,e) as a function of e has the Taylor expansion: 

oo 

uj^'\t,e)=J°\t) + Y,e''J''''\t), te[0,T], (5.12) 

n=l 

where 

The function a;*^^'"^^ satisfies the equation: 

^{s,n) ^ v/(a;(°^) o + /i^^'") , i G [0, T] , (5.14) 

where 



nl \ 



1 / a" 



/(a;(^)(t,e)) + e^(a;(^)(t,e)) 



- V/(a;(°))oa;(^'") 

e=0 
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for example 

^(.,2) ^ l|vV(u;W)oa;(^'i)|ou;(^'i) 
+V^(u;(°))oa;(^'i) , 
= |vVM°))ocu(^'i)|oa;(^'2) 

+V^(a;(°))oa;(^'2) 

Next we are going to derive an expression for the time derivative of the functions A[^'"^ 
= (V^t/ o / , a;(^'")) + (VC/ , V/ o a;^*'")) + (Vt/ , /i^^'")) . (5.15) 

Notice that 

(VC/,/) = 0, (5.16) 

from the invariance of U under the e = flow (4.1); then the variation of the equation 
(5.16) leads to: 

(VV oSu, f) + (VC/ ,Vfo5u;)^0, 
and since V^C/ is a symmetric operator, we have 

{V^U o f , 600) + {VU ,Vf oSoo) = 0. (5.17) 

If we substitute 5uj in (5.17) by a;^*'"\ we see that equation (5.15) is reduced to 

AJ^''^^ = {VU , , t G [0, T] . (5.18) 

Without rigorous justifications and taking the limit T +00 in (5.18), we have 

A[;'")(0) = - /"(Vt/(cu(°)) , dt . (5.19) 



Similarly, 

A[7'"^(0) = /° (VC/(a;(°)) , /i^"'")) , (5.20) 

J— 00 

where has the same expression as with s replaced by u everywhere, 

00 

A^{t,e) = E6'^A[7''^)(t) , te(-oo,0], (5.21) 

n=l 
1 Qn 
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= ^^c^m , g-ry'^Kt,o)) , (5.22) 

oo 

J''\t,e) = u;(°Hi) + 5]e"a;("'")(^) , ie(-oo,0], (5.23) 

n=l 

1 8^ 

= ^^u;(")(i,0) , (5.24) 
^{u,n) ^ V/(u;(°^) o + , i e (-00, 0] . (5.25) 



From the expression (5.9) of the signed distance du, we have 

oo 

du = 

n=l 
oo 

n=l 

where M^"^ is called the n-th order Melnikov function and M^^ is the usual Melnikov 
function. 

Notice that cu(")(0,e) and cj(*)(0,e) are chosen to have the same coordinates except the 
VU{0) and VV'(O) directions. From the expressions (5.10) (5.11) (5.12) (5.13) (5.21) (5.22) 
(5.23) and (5.24), if M^'^ = 0, j = 1,2, ■ ■ ■ ,n - 1; Y = U,V; then 

J''^\0) = u;("'^')(0) , J = 1, 2, • • • , n - 1 . 

Notice that both a;^*''-') and o;*^"'-'^ satisfy the same form of equations (5.14) and (5.25); thus 
(^(sJ) and o;^"'.'^ together represent one solution on the entire interval t e (— oo, oo), and we 
can drop 's' and 'u\ and use the simple notation u^^^ for j' = 1, 2, • • • , n — 1. 

+ /i(^) , t e (-00, oo) , (5.26) 



A^-'"^(0)-A^+'"^(0) 



j = 1,2, •••,n - 1 , 



where 



f{u^^\t,e)) + eg{u;^^^{t,e)) 



-V/(u;(°))oa;(^) , 



e=0 



(5.27) 



l,2,---,n , 



for example 



^|vV(c^(°))oa;«|oa;« 



,(1) 
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+1 1 1 V V M° V ^^'^ } ° ^^'^ } ° ^^'^ 

+ i|v25(a;(°))ou;«|oa;« . 

Notice that both e) and Cij^^\t) approach the same point —u;* in backward time, we 

have 

^(«.0(_oo) = 0, Vi = l,2,---. (5.28) 
Prom the expressions (5.19) and (5.20), we have 

J—oc 

similarly for My\ 

Summciry. The signed distances dy {Y = U,V) defined in (5.7) and (5.8) have the 
representations: 

oo 

dy^Y. ^"^y ^ ' {Y = U,V) . (5.29) 

n=l 

IfM'^^ = 0, J = 1, 2, • • • , n - 1; y = [/, V ; then 

, , /'OO 

mP = / (Vr(u;(°)) , , F = [/, ; 

J — oo 

where /i*^"-* ^'wen m (5.27) and uj^^^ (j = 1,2, ■ ■ ■ ,n — 1) solves the linear equation (5.26) 
under the boundary condition (5.28). 

5.3 Numerical Evaluations of the Second and Third Order Mel- 
nikov Functions 

5.3.1 Numerical Evaluations of the Second Order Melnikov Functions 

In the formulae of the second order Melnikov functions [(5.29), for n = 2], since the functions 
U and V only depend on (x'j (1 < j < 4), only the following four components of /i*-^-* enter 
the evaluation: 









(5.30) 








(5.31) 








(5.32) 






■Ao[a;i°)4na;Wa;f]. 


(5.33) 
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Thus we need to know the functions Oi)^^\ Uj^^ (0 < j < 5). Specifically these functions 

satisfy a set of decoupled systems of linear equations. The five functions ou^\ ou^^^ (1 < 
£ < 4) satisfy a self-contained system of linear equations with variable coefficients: 



/ 



UJ 



(1) 



^3 



I 



A 



(1) 



(1) 



(1) 



(5.34) 



where 



A 



B 
















(0) 



Ao4'^u;^o 





(0) 







,(0), ,(0) . , ,(0), ,(0) 









-A2U)\ 



(0) 



,(0) 



(0) 



2t(;2°'' — A\ijj, 

A,uf 




(0) 



(0) 

•1 



The functions (^q"*^^ and Jf^ satisfy the hnear equations: 



- ujf\A,J^-A,u^f\^A,.^f^f. 



(5.35) 
(5.36) 



The two functions uj^^\ and uo^'' on the right hand sides of (5.35) and (5.36) are obtained 
from solving the following two self-contained systems of linear equations: 



(1) 



d_ 

dc 



UJ 



CJ_3 



c 



. ,(1) 

a;_2 



(5.37) 
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d 
dC 



where ^ jojft^ ^ = ^ cosh{KTt + tq), 



C 







^6 





















(5.38) 



V 

















A 3 


























r 

^ 82 










23 Q 

50 " 


59 

12s 


- 




^ 

^ 82 





83 
170 







59 
122 







" -A„c.r - 




26 '^O 


































-A7 










Ae 


-^8 




Ay 





-A, 







A, 








-c, 



^50;^' 






(0) 



11,, (0) 
~26'^5 







The characteristic equation for the matrix C: 
"23 X 39 



+ 



59 ^39 83 
~ ^ 170 



.50 X 82 ' 122 V82 
has only imaginary eigenvalues, 

Ai,2=±i 0.773, 



+ 



23 X 39 X 59 X 83 
50 X 82 X 122 X 170 



0, 



^3,4- 



±i 0.295. 
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Let Vj (1 < J < 4) be the corresponding eigenvectors, v — columns {vi,V2,Vs,V4}, 
w~ — a;i^2) <^-3) '^-i)'^) '^'^ — (<^6^\ <^7^\ <^8^\ '^g^'*)'^; then 

Jsign(/t)oo I J 

v-^V{i) , (5.39) 

Jsign(K)oo L J 

^-'^(0 ■ (5.40) 

The solutions w"^ given above satisfy the boundary condition w'^lt — — oo) = 0. In fact, 
w"^ vanish at both positive and negative infinities, w'^{t — oo) = 0. 

The system of coupled linear equations (5.34), (5.37) and (5.38) were integrated nu- 
merically by Thomas Witelski using a fourth order Rungc-Kutta scheme. The resulting 
numerical solutions were then used to evaluate the integrals for the second order Melnikov 
function Mjj^ and My^ using the trapezoidal rule. For calculations with sufficiently small 
step-sizes, the second order Melnikov functions were found to be identically zero to machine 
precision, independent of the value of the parameters 6o and tq. Evaluation of the partial 
integrals for — oo < t < t strongly suggest that the integrands of M^^ and My^ are odd 
functions. 



5.3.2 Numerical Evaluations of the Third Order Melnikov Functions 

In the formulae of the third order Melnikov functions [(5.29), for n = 3], since the functions 
U and V only depend on o^j (1 < j < 4), only the following four components of h^^^ enter 
the evaluation: 



(3) 



Uq — 



A2UJ, 



(1) 



U)2 



+ 



A^Lof^ - A2UJ2 



(1) 



(2) 



+ 



(1) 



UJ. 



(2) 



L02 



AlU, 



(1) 



(1) 



(2) 



(2) 



UJ 



(2) 



+ 



(2) 



Thus we need to know the functions uj'^\ a;j ' {0 < j < 5). Specifically these functions 
satisfy a set of decoupled systems of linear equations. The five functions 
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£ < 4) satisfy a self-contained system of linear equations with variable coefficients: 



(2) 



^3 



(2) 



(2) 
1 

(2) 



^3 



(2) 



where A is given in (5.34), and 



where hf\ h^2 \ hf^ ^^'^ h^i' ^-^^ given in (5.30), (5.31), (5.32) and (5.33), and 



,(2) 



— A_4__3a;_4a;_3 + A._^_2'^-z^-2 + ^-2,-i<^-2'^-i 

+^-l,0'^-l'^0 +^0,ll'^0 "^l +'^0 ) 



-Ai^2^^1 1^2 + ^2,3^^2 ^3 



^3,41^3 1^4 



+A4,5(o.f4^^ + 4^^-f)+A5^ac.f4^) 
_l4 , _l 4 ,,(1), ,(y_L4 ,,(1), ,(1) 



(1) 



This expression is obtained from the fact that 0;^ 
functions oJq'^ and ojf^ satisfy the linear equations: 



for / < -5 and I > 10. The two 



,(2) 



(2) 



UJ 



(2) 



/,(2) 



where 



Ci^4 -|- 



A4a;4 



(0) 



(2) 



(0) 



UJ. 



(2) , /,{2) 



<^6 + "5 ) 



(5.41) 
(5.42) 
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The two functions u^^l and on the right hand sides of (5.41) and (5.42) are obtained 
from solving the following two self-contained systems of linear equations: 



d_ 



d_ 



where (, C and S are given in (5.37) and (5.38), 



p(2) _ 



(2) 
-1 

(2) 



r (2) - 




r (2) 1 




('2^ 


= C 


tc;_2 












r r 1 




^6 




^8 


= 8 


^8 








<^9 





^(2) 



/,(2) 



(2) 
7 

(2) 



(5.43) 



(5.44) 



h 
,(1) 



(2) 



(1),„(1) 
P 



Ao(4«)4na;«4°^), 



A-3UJ. 
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In these formulas, u^l and a;io'* appeared, and they can be obtained from. 



(1) 
10 



A I Ml 
— A_4a;^ ''a;_4 , 



a; 



(1) 



cui^i'' , cji^T , iJ3x , c<j^] and 



Let 

representations as for w~ and w'^ ^ 



w-{C) 



t diag I 

J sign(/t)oo I 



(4'\4'\4'\4'y; then we have the similar 



gAi(c-0 e^3(C-0 gA4(C-0 
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v-'V^'\0 , (5.45) 
C 



Jsign(K)oo L J 

v-'j^^'\Od^- (5.46) 



The solutions given above satisfy the boundary condition w'^lt — — oo) = 0. In fact, 
w"^ vanish at both positive and negative infinities, = oo) = 0. 
The above systems are solved under the boundary condition, 

6uP(t = -oo) = , yiez. 

(2) 

Moreover, direct verification shows that LOf = for / < —10 and / > 15. For this case, 
the numerics of Thomas Witelski could not give a firm conclusion, and we are continuing 
investigating the numerics. This numerics is a part of a future project on numerical inves- 
tigation of the degeneracy v.s. nondegeneracy of the hyperbolic fohations of the 2D Euler 
equation. 



6 Conclusion and Discussion 

In this Part II of our study, we see an explicit representation for the degeneracy of the 
hyperbolic structures for a Galerkin truncation. We use higher order Melnikov functions 
to study the robustness of such structures for the so-called dashed-hne model. We find 
that both the first and the second order Melnikov functions are identically zero, which 
indicates that such structures arc relatively robust. The study in this paper serves as a 
clue in searching for homoclinic structures for 2D Euler equation. The recent breakthrough 
result [24] of mine on the existence of a Lax pair for 2D Euler equation makes it plausible 
for the existence of homoclinic structures. 
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